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Abstract: We study Coleman-de Luccia tunneling in some detail. We show that, for a 
single scalar field potential with a true and a false vacuum, there are four types of tunnel- 
ing, depending on the properties of the potential. A general tunneling process involves a 
combination of thermal (Gibbons-Hawking temperature) fluctuation part way up the bar- 
rier followed by quantum tunneling. The thin-wall approximation is a special limit of the 
■ case (of only quantum tunneling) where inside the nucleation bubble is the true vacuum 

while the outside reaches the false vacuum. Hawking-Moss tunneling is the (only thermal 
. fluctuation) limit of the case where the inside of the bubble does not reach the true vacuum 

at the moment of its creation, and the outside is cut off by the de Sitter horizon before 

it reaches the false vacuum. We estimate the tunneling rate for this case and find that 
00 

the corrections to the Hawking-Moss formula can be large. In all cases, we see that the 
bounce of the Euclidean action decreases rapidly as the vacuum energy density increases, 
signaling that the tunneling is not exponentially suppressed. This phenomenon may be 
interpreted as a finite temperature effect due to the Gibbons-Hawking temperature of the 
de Sitter space. As an application, we discuss the implication of this tunneling property to 
the cosmic landscape. 
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1. Introduction 

Tunneling in the presence of gravity was first beautifully analyzed by Coleman and de 
Luccia (CDL) This problem emerges in the study of the early universe in relation to 
the inflationary epoch, and more recently, in relation to the cosmic landscape as suggested 
by superstring theory. However, there are properties of the CDL tunneling that remain to 
be better understood. 

Consider a simplified theory of a single scalar field with potential V((f>) and a canonical 
kinetic term in the presence of gravity. The potential V((p) has a false vacuum V+ at 0+ and 
a true vacuum V- at 0_, so the energy density difference between them is e = V+ — V- > 0. 
The potential barrier between these two local minima has height Vr, as shown in Figure |l|. 
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Figure 1: Potential V (</>). 



In the semiclassical approximation, the tunneling rate per unit volume is given by 



T ~ Ae 



-B 



B = S E {cj>) - s E (h 



(1.1) 



ri.2 



where Se(<P) is the Euclidean action for the bounce solution and Se(((>+) is the Euclidean 
action evaluated at the false vacuum at </>+. In the thin- wall approximation (e — ► 0), when 
the vacuum energy density is negligible compared to the Planck scale, B recovers the result 
in the absence of gravity ||], as expected. On the other hand, as we increase V + > V_, 
while still in the thin- wall approximation 01, B takes a different form, 



B 



27vr 2 r 4 2vr 2 r 



2e 3 



H 3 



(1.3) 



where r is the domain wall tension and the Hubble constant is given by H 2 = 8irG]\rV+/3 = 
V+/3Af 2 . Note that the first form (when V+ is small) is independent of V+ while the second 
form is independent of e = V+ — V_. We see that, for fixed r, B decreases rapidly as H 
increases. It seems that we can easily have a situation where B <C 1, in which case, the 
tunneling is not exponentially suppressed at all (and the evaluation of the prefactor A in 
Ql.l| ) as well as sub-dominant contributions to T becomes important). 

A similar phenomenon seems to be happening in the Hawking-Moss (HM) tunneling 



B HM = 2^ l Mi 



1 



l 



8irAV + 
3# 4 



(1.4) 



where AV + = Vt — V+- Note that if we move V(4>) up without changing its shape, that is, 
keeping e, r and AV+ fixed, B decreases like Q r H~ 4 , depending on which formula (i.e., 
thin wall or HM) is applicable. In this paper, we like to see if and under what conditions 
this phenomenon is real. More precisely, we examine the CDL tunneling more generally to 
see when the thin-wall approximation or the HM formula ( |1.4[ ) is valid. 
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CDL tunneling concerns quantum fluctuations in the false vacuum of a nucleation 
bubble which subsequently grows classically. The following picture emerges. Depending on 
the properties of V((f>), there are four regions in the parameter space, yielding four different 
situations how the nucleation bubble is created : 

(I) The center of the nucleation bubble is in the true vacuum V- while the outside of 
the bubble reaches the false vacuum V+. The thin- wall approximation is a special limit 
(e — > 0) in this case. 

(II) The outside of the bubble reaches the false vacuum V + but the inside of the 
nucleation bubble never reaches the true vacuum V— in the Euclidean solution. After the 
creation of the bubble, its inside falls towards the true vacuum V- as the bubble grows. 

(III) The inside of the bubble reaches V—, but the outside of the bubble does not reach 
the false vacuum in the Euclidean action, due to the presence of the de Sitter horizon. 

(IV) Not only does the inside of the bubble not reach V- , the outside of the bubble 
does not reach the false vacuum in the Euclidean action, due to the presence of the de 
Sitter horizon. HM tunneling is a special limit in this case. 

In the absence of gravity, regions faraway from the bubble are by definition in the false 
vacuum, so we have only case (I) and (II). In the presence of gravity, the de Sitter horizon 
cut off the scalar field 4> as it tries to reach the false vacuum at large distance from the nu- 
cleation bubble pi. That is why case (III) and (IV) are possible. This phenomenon of cases 
(III) and (IV) may be interpreted as a Gibbons-Hawking (GH) temperature Tjj = H/2tt ef- 
fect. This point was recently made by Brown and Weinberg ||, who show that tunneling in 
the presence of gravity may be interpreted as a combination of thermal (Gibbons-Hawking 
temperature) fluctuation plus quantum tunneling. We agree with this picture. That is, (f) 
thermally fluctuates from the false vacuum at part way up the potential (to </>/+) before 
quantum tunneling to the other side of the barrier. In the Hawking-Moss limit, it fluctuates 
all the way to the top of the barrier at Vt- Here, our explicit calculations (in particular, 
the determination of </>/+) allow us to estimate in some detail, for a given potential, the 
individual contribution of the thermal fluctuation versus that of the quantum tunneling in 
the tunneling rate. We also evaluate the back-reaction effect on the background geometry 
and show that it can be very important. 

There are a few other properties that are worth mentioning here. 

• The thin- wall approximation is valid only when the tunneling belongs to case (I). In 
addition to the case (I) conditions (to be specified), we see that Mp > Vr > AV+ S> e. 
As we increase Vt with a generic fixed shape of V((f>) (with e > 0), the tunneling goes over 
to case (IV) via case (II) and/or case (III). 

• For large Vt and e > 0, case (IV) is generic, where HM tunneling is a special limit. 
Here we estimate the accuracy of the HM formula (1.4) and its correction. We see that the 
HM formula for B is an over-estimate, by as much as a factor of three. That is, the actual 
decay is a combination of thermal tunneling and quantum tunneling, whose rate can be 
much faster. 

• Another interesting point is that tunneling in the absence of gravity is always down- 
wards, while gravitational effects allow tunneling upwards. 
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Figure 2: The log-log plot of the tunneling exponent B as a function of the Hubble parameter H 
for a fixed potential V(cj>) except for its overall height as measured by H. (G = 1). The three cases 
are separated by the two solid vertical lines. As H ~ \/Vr increases, the tunneling behavior goes 
from case (I) to case (IV) via case (II); that is, it goes from the CDL thin- wall approximation (to 
the left in case (I), with slope ~ —3 as given by Eq.(1.3)), to the HM tunneling (to the right of 
case (IV), with slope ~ —4 as given by Eq.(1.4)). For this particular potential, the slope case (IV) 
(~ —3.74) never quite reaches the HM value due to corrections. The specific potential used here is 
described in Sec. 6. Case (III) does not appear for this potential. 



• With an overall picture, we see that the fast drop off of B as H increases is real. This 
is illustrated in Figure |||, where case (III) does not appear due to the particular choice of 
the potential used. 

• Strictly speaking, for a "general" smooth potential with varying parameters, the 
above four cases may reduce to case (IV) only. However, within case (IV), it contains 
regions that resemble the various cases discussed above. The triangular potential has four 
different cases so it actually demonstrates the change from thin wall to HM case more 
dramatically (but does not change the physics qualitatively.) 

• As B decreases, multi-nucleation bubbles will form and the phase transition involves 
bubble collisions. Once B is small (say B < 1), Eq. ( |1 . 1| ) is no longer valid. The prefactor 
and sub- leading terms will become important. All we can say is that tunneling is no longer 
exponentially suppressed. The rate of transition depends on the details. We do expect 
that 0(4) nucleation bubbles no longer dominate. That is, bubbles of other shapes as well 
as bubble collisions become important. The transition can be complicated and is at a fast 
time scale. Something analogous to spinodal decomposition may happen. 

In terms of pure thermal (Gibbons-Hawking (GH) temperature) tunneling the sup- 
pression of the tunneling is interpreted as due to the Boltzmann factor, where the bubble 
is a 3-sphere and the inverse temperature is treated as an imaginary time with period 
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1/Ttf. In case (I) and (II), we see that the quantum tunneling is dominant (i.e., has a 
smaller value) due to the 0(4) symmetry that is lacking in the finite temperature effect. 
So strictly speaking, the enhancement of the tunneling rate due to a large vacuum energy 
density is a pure gravitational effect. In de Sitter-like vacua, we may interpret this as a 
GH temperature effect in the presence of an enhanced (here 0(A)) symmetry. In case (III) 
and (IV), the GH temperature starts to play a more prominent role. 

The same GH temperature will also contribute to the finite temperature effect on the 
potential V(4>) — * V(4>,Th) [jj], ||. This contribution is perturbative in the couplings while 
its effect in tunneling is non-perturbative. In the finite temperature formalism, this leads to 
a term of the form T^cft 2 into the finite temperature potential. In the gravity perspective, 
this is simply a coupling of the form R(p 2 , where R is the Ricci scalar. One an also interpret 
this as a finite volume effect, due to Gauss's law, H 2 (j) 2 , where H~ l is the horizon size. 
The overall picture is self-consistent and clear. 

We shall comment on the impact of the GH temperature on the cosmic landscape, both 
on the shape of the effective potential of the landscape and on the tunneling rate. Because 
of fast tunneling when the wavefunction of the universe is high up in the landscape, it is 
likely that the universe is quite mobile there. 

For a general potential V((p), the coupled equations of <?!>(£) and the cosmic scale factor 
r(£) are too complicated to solve even in the absence of the gravity except in the thin- wall 
approximation. However, for some special potentials, say, the triangular potential V((j)) 
(Figure ||), the bounce B can be obtained analytically |J. So it is natural to consider the 
tunneling rates for such a potential in the presence of gravity. Unlike Ref.Q, we cannot get 
an exact analytic formula for B in the presence of gravity, but the resulting analytic study 
gives a very good approximation and does simplify enough for us to see the overall picture. 
Here the absolute height of the potential, which corresponds to the vacuum energy density, 
is important for tunneling with gravity. To be specific, we shall adopt this triangular 
potential in this paper. 

The rest of the paper is organized as follows. In Sec. 2, we briefly review the overall 
framework of tunneling in de Sitter space. This framework is the CDL tunneling formalism. 
In Sec. 3, we go back to tunneling in the absence of gravity. In particular, we review 
the special case of tunneling in a triangular potential studied by Duncan and Jensen ||. 
Here we emphasize the formulation that is suitable in the extension of their analysis to 
include gravity. Sec. 4 presents the setup for tunneling in de Sitter space, again using the 
triangular potential. As 4> varies, the tunneling is not happening in pure de Sitter space and 
the back-reacrtion is estimated. Sec. 5 presents the main result of this paper. The above 
four cases and their conditions are discussed. The meaning and implications of the results 
are discussed in Sec. 6. We then discuss thermal tunneling in Sec. 7. Here we are referring 
to the treatment of the de Sitter horizon effect as a Gibbons-Hawking temperature effect. 
In Sec. 8, we point out that the Gibbons-Hawking temperature should also modify the 
potential via finite temperature effect on potentials. We then consider the implication of 
the Gibbons-Hawking temperature on the cosmic landscape. Sec. 9 contains the summary 
and some remarks. 
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2. Coleman-de Luccia Tunneling 



Let us consider the theory of a single scalar field 
of gravity, given by 



with a potential V(4>) in the presence 



S 



-9 



(2.1) 



where k = 8vrG = 1/M 2 . 

The potential V(4>) has a false vacuum V+ at + and a true vacuum V- at <j>— There is 
a potential barrier between these two local minima, as shown in Figure |l[ Let the height of 
the barrier at ^ be 1^. We have chosen 0_ > fix > 4>+- Let the energy density difference 
between the false and the true vacua be e = V+ — V- > 0. 

In the presence of gravity, the CDL tunneling rate per unit volume is given by Eq. 
(1.1. l~2l) , in term of the coefficient B = Se{4>) — Se{4>+)i where Se(4>) is the Euclidean 
action for the bounce solution and Se(4>+) i s the Euclidean action evaluated at the false 
vacuum. Since a Euclidean solution with a 0(4) symmetry has an etremum action, B is in 
general dominated by the "bounce" solution with 0(4) symmetry. This solution has the 
Euclidean metric, 



ds 2 



de + m 2 dnl 



(2.2) 



where dQ 2 is the metric of a unit 3-sphere. The Euclidean equation for the "bounce" 
solution is determined by the minimum value of the Euclidean action, 

3r' 

r 

2 ( 



+ 



„'2 



1 + 



3M 2 



f/2 



dV 



V 



(2.3) 
(2.4) 



where the prime is derivative with respect to £. Here the Einstein equation yields one 
non-trivial equation, since the equation of motion for r(£) follows from Eqs. ( |2.3 |Q ). We 
can choose r(0) = 0. Using Eq.(2.4) to simplify the Euclidean action, one obtains 



S E = 4vr 2 / d£ 



r 6 V 



(2.5) 



The CDL instanton is a unique solution with the topology of a four sphere S 4 in 



Euclidean space |10|. The geometry after bubble nucleation is described by the analytic 
continuation of the CDL instanton to Lorentzian signature. The radial coordinate £ is 
continued to £ = it and the metric in Lorentzian frame is 



-ds 2 = -dt 2 + r 2 (it)dQ 2 H3 , 



(2.6) 



here the metric is multiplied by an overall minus sign and dfl H 3 is the element of length for 
a unit hyperboloid with timelike normal. The metric within the spherical bubble describes 
a spatially open Friedmann-Robertson- Walker Universe. 

As we shall see, a qualitative picture emerges: depending on the properties of the 
potential, there are four different cases how the nucleation bubble is created : 
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(I) The center of the nucleation bubble reaches the true vacuum V- while the outside 
of the bubble reaches the false vacuum V+. The thin- wall approximation is a special limit 
(e — > 0) in this case. 

(II) The outside of the bubble reaches the false vacuum V+ but the inside of the 
nucleation bubble never reaches the true vacuum V- in the Euclidean solution. After the 
creation of the bubble, its inside will fall towards the true vacuum V- as the bubble grows. 

(III) The inside of the bubble reaches V-, but the outside of the bubble does not reach 
the false vacuum V+ in the Euclidean action, due to the presence of the de Sitter horizon. 
This case and the next case never happen in the absence of gravity, since we start with 
the false vacuum everywhere. Far from the bubble, eft approaches the false vacuum by 
definition. However, when the vacuum energy is not negligible, there is a de Sitter horizon 
so it is possible that eft hits the horizon before it reaches the false vacuum. 

(IV) Not only does the inside of the bubble not reach V- , the outside of the bubble 
does not reach the false vacuum in the Euclidean action. The HM tunneling is the limit in 
this case. Here we are able to estimate the accuracy of the Hawking-Moss formula and its 
correction, which can be substantial. 



For a general potential V{4>), the coupled equations (2.3, 2.4 ) are complicated to solve. 
In the absence of the gravity, for some simple potentials, say, the triangular potential V(4>) 
(Figure |3|) , the bounce B can be obtained analytically || . Then it is natural to consider 
the tunneling rates for such potentials in the presence of the gravity. Unlike Ref.||, we 
cannot get an analytical formula for B in the presence of gravity, but the analysis does 
simplify enough for us to see the overall picture. Here the absolute height of the potential, 
which corresponds to the vacuum energy density, is important for tunneling with gravity. 
To be specific, we shall adopt this triangular potential in this paper. 

The triangular potential may be parametrized in the following way. Let the height 
of the barrier of V(4>) at <jyp be Vt, which also provides a measure of the overall vacuum 
energy density. Let AV± = (Vt — V±). Both Acp + = 4>t — 4>+ and = 0_ — 4>t are 
defined to be positive so the slopes (gradients) ±A± of V(4>) are given by 

A<p± 

Note that, in general, A + ^ A_. Since the e = AV_ — AV + > 0, we have 



AV+ _ 1 

v V AV - V 1 + t/ Av - 

Acf> + , A0_ and A. Here, to simplify the solution, we set A^ = Xr = A. 



3. Tunneling Without Gravity 

When gravity is negligible, Eq. (|2,4j ) reduces to r' = 1. This happens if we set k = 0. 
Alternatively, this is a very good approximation when the overall height of the potential is 
much smaller than the Planck scale. 
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Figure 3: Triangular potential V(<f>). The false vacuum is at V+ = V(<fi+) and the true vacuum is 
at V- = V ((/>-), with the top of the barrier at Vr — V{4>t)- Here A+ and — A_ are the gradients. 



With r(£ = 0) = 0, we have r = £, so tunneling reduces to the simple case without 
gravity. The solution must satisfy the boundary condition, 



0'(O) = 



(3.1) 



which is necessary for <f> to make sense at the center of the nucleation bubble. At large 
radius, we expect 

<^0+, £^oo (3.2) 

As we shall see, this condition may be modified when gravity is important. 

Let us first review the solution to the triangular potential without gravity, the case 
worked out by Duncan and Jenson M . For <p + <</><</>- , Eq. (EO) becomes 



(3.3) 



The general solution is 



m = a+ b/e ± A±e 2 /s 



(3.4) 



where the constant a and b are determined by the boundary conditions and by matching 
the field values and their derivatives at the top of the barrier, which occurs at some radius 
£t to be determined. 

It is easy to see that cj) will reach its false vacuum value at finite radius (to be 
determined) and then stay there. The above boundary condition (|3.2| ) is replaced by 



(3.5) 



while there are two possibilities to satisfy the condition (|3.1|), The first case, namely case 
(I), is when 4> stays close to until at radius ^_ when <p starts to decrease under Eq. (EO) . 
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The boundary conditions in this case are 

j m = <i>- o < { < {_ 

1 *'(«-) = o (3 ' 6) 

The other possibility, namely case (II), happens when, at the time of creation, the inside 
of the bubble never reaches the true vacuum. In this case, the boundary conditions are 

/ 0(0) = <Po r „ 7 N 

\ 0'(o) = o (3 - 7) 

where the initial value 4>o is to be determined. This case only happens if (fio is to the right 
side of the barrier (i.e., near </>_) and cfto < </>_. Also, we expect V + > V{4>q) > V-. 

3.1 Case (II) 

Let us consider case (II) first. On the right and left sides of the barrier, we have 

f mo = 4>o-^e o < e < 
\ mo = <p+ + - e + ? tT<t<u 

Matching the derivatives of the two solutions Q3.8] ) at 

£ = (1 + c)& (3.9) 
where c = A_/A + . Matching the field values at £t yields 



(3.8) 



(3.10) 



Now that the unknowns £t, and (fio are determined in terms of properties of the potential 
y (</>), we simply insert the solution ( |3.8| ) into the Euclidean action and integrate it from 
£ = to ^ = £+. The bounce (1.2) is given by 

B = ^ J^ C ^ (3.11) 
3 (^TT^-1) 4 AP+ v ; 

Let AV = Vt-V = V t - V(#o) so, using Eq.fl3T0D, 

Ay+ ^ Ay + _ (vT+c- 1)" 



Ay_ < a£ = c2 (3 - 12) 



which implies that tunneling is always downward, and A V+ / A Vq < 1 /4. 

Recall that case (II) holds only if Vq > V- , or <j)Q < (/>_ , which translates to the above 
condition (3.12). or equivalently, /3 > 1, where 



Ay_ A0--A</> + 

2A0- (3 - 13) 

In general, /? is semi-positive. If j3 = 1, then 0o = 0- an d one can rewrite the bounce 
(3.11) in the following form 

B= *W--**)' (3 . 14 , 
3 AV+ V ; 
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3.2 Case (I) 

If j3 < 1, the inside of the bubble reaches the true vacuum V_, that is, </>(0) =</>—. At 
certain radius begins to decrease until it reaches 4> + at radius £+. The region for 
£ > £_l stays at the false vacuum. Now, the solution </>(£) contains four pieces, 



m = { 



MO 
h(0 



+ 



8£ 2 



(3.15) 



where <j)- > </>i(£) > (pT and 0+ < </>2(£) — Now we have three unknowns : ^t, £- and 
£+. Matching the derivatives as well as the field values at £t yields 



A0_ 



2\2 



(3.16) 



(e 2 



T 



so ^T) £- and £+ can be solved in terms of the properties of V{4>). Once again, one can 
insert these solutions of the parameters into the 4> solution ( 3.15| ) and evaluate the bounce. 
We find that it is convenient to introduce the tension to be 



T+ + T_ 



r<t>T o 
r± = ± / d^2[V{<i)) - V±] = - V / 2A^A<A± 

U± 6 



After some algebra, B can be written in terms of these values 



B 



27n 2 
32? 



[r+(l + 1/t/) + r_(l + r?)] 3 [r + (3/r ? - 1) + r_(3r/ - 1)] 



(3.17) 



(3.18) 



(3.19) 



1, .B in ( p. 19 ) reduces 



where e = V + - V- and rf = AV + /AV- < 1 (fPp . In the limit /? 
to B in ( |3.14| ), that is, it agrees with B in case (II). 

Here t± in Eq. ( |3.18| ) should be treated as a convenient definition. The value of r as 
defined should be very close to the actual domain wall tension. For more general potentials, 
the coefficient 2/3 in Eq. ( |3.18| ) will have order one variations. Note that e/ AV- = 1 — rf . 
In the thin- wall approximation, rj — > 1 as e/AV- becomes small, and B reduces to the 
usual thin- wall formula for B, 

.-,--2 

r 4 (3.20) 



B 



277T 2 

2e 3 ' 



So we see that the thin- wall approximation is inside the region of case (I). Comparing the 
small e/AV- case to that in case (I) Eq.(3.12), we see that there is a sizable parameter 
region where the thin- wall approximation is not valid. For case (II), since the inside 
of the bubble never reaches the true vacuum value </>_, the wall tension does not carry 
much physical significance here. The non thin- wall approximation has also been studied in 
Ref.|ll. 
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In summary, the regions (I) and (II) are divided by the value of ( 3,13| ): case (I) if 



< < 1 and case (II) (<p(0) = 4>o < 4>-) if > 1. For = 1, only the center of the 
bubble reaches the true vacuum. 

In the absence of gravity, the overall height of the potential is not important, so the 
triangular potential is parameterized by a set of four parameters, namely (AV±, AA±). 
This set can be replaced by an equivalent set (AV±, A</> + , c). The bounce B in Eq. Q3.11| ) 
is expressed in terms of this set, where we see that AV- does not enter, since the inside 
of the bubble never reaches V-. In Eq.( |3.19| ), we express B in terms of an equivalent set 
( T ±> e > r i)i which is more convenient when we want to take the thin- wall limit. Here, we 
also see how B behaves as we move away from the triangular potential: the leading order 
correction is automatically incorporated into the two tension components r± by varying 
appropriately the 2/3 coefficient in Eq.( |3.18D . 

If A_ = A+, we see that simplifies somewhat, 

( , 21) 

where < rj < 1 is given in Eq.plf); so the dividing point is at r] c = (y/2 - l) 2 = 0.17. 
Often times in our analysis below, we shall restrict ourselves to this special case. 

4. Turning on Gravity 

Let us now turn on gravity Q, so the triangular potential has five parameters, namely 
(AV±, AA±, Vt)- As noted above, at times, it may be convenient to choose an alternative 
but equivalent set of parameters. 

To emphasize the effect of gravity, let us consider the situation where 

V T > AV- = V T -V- (4.1) 

where Vr/Mp C 1 so gravitational effects can be important for the tunneling while the 
semi-classical approximation is still valid. 
The equation fl2.3| ) now reads, 



where L = <j>' 2 /2 — V + Vr is the (shifted) Euclidean energy. By fl4.1|) , for large Vt, we use 
the approximation 

r' 2 ~ 1 - -^V T = 1 - H 2 r 2 . (4.3) 

3M p 2 v ' 

where 

H 2 = Vr/3M 2 (4.4) 
is the Hubble constant of the de Sitter space with vacuum energy Vt- So we have 

r(f) = H- 1 sm(H^) (4.5) 
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where r(0) = r(ir/H) = 0. That is, the range of £ is bounded, < £ < tt/H. 

It is easy to check that this approximation is self-consistent, i.e., ^ Vt- For a 

tunneling solution, 0(£) G so 

L > 0. (4.6) 

Using Eq.Q, 

L' = Q<// 2 - V(</>)\ = -3^0' 2 ~ -3Fcot(FO0' 2 (4.7) 

When £ G [0,7r/2iT), L' < and if £ G (ir/2H,ir/H], V > 0. Therefore maxL = 
max{L(0),L(7r/#)} < AVI. So \L\ < Vr, and Eq.(|4~5|) is valid to the leading order 
in AV-/Vt, For later purposes, we note that 

24tt 2 M^ 

S E (<P + ) = ^ ( 4 -8) 



The equations (2^3, 2A) are now decoupled, and 



cP" + 3Hcot(H£)cP' = (4.9) 
d(p 

and can be solved with appropriate boundary conditions. This we shall do in the next 
section. 

The bounce solution 0(£) will modify the geometric background used, i.e., r(£) in 
( |4.5| ). As we shall see, it is important to include this back-reaction. Inserting 0(£) back 
into Eq.( [4.2D , we get a first-order differential equation of the leading-order correction of r, 
namely <5r(£), 

cos(H£)6r' = -Hsin(H£)6r + sin 2 (#£). (4.10) 
The formal solution for 5r is, 



5r(0 = cos(H0 ■ / ? dr/tan 2 ^)^, (4.11) 
Jo 2 Vt 



where the superficial singularity at rj = ir/(2H) can be regularized as 

l( v )-l(^) 



6r{£) = cos(H£) ■ dr) tan 2 (Hi]) 
Jo 



2V T 



2Vt 2H 1 ' 1 ' J 
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Notice that by ( ^4.7|) , L'(ir/2H) = and the integral is well defined. Now we can consider 
the Euclidean action for this solution, to the leading order of AV-/Vt, 



S e = 4tt 



+ 



d£( (r + 5r) 6 V(£) - 3(r + 5r)M^ 
At: 2 j 1 d£ (r 3 V T - SrM$) 
4tt 2 ' d^Sr(3r 2 V T - 3M 2 ) + r 3 (V(£) - V T )\ 



+ O 



AV 2 \ 



(4.13) 
(4.14) 
(4.15) 



We separate the result into three parts. The first part ( 4.13[) is — 247r 2 Mp/Vr, the de-Sitter 
space Euclidean action. The second part (|4. 14|) is suppressed by a small factor AV-/Vr 
comparing with the first part, because by ( |4.12| ) the magnitude of 5r is of H~ l AV_/Vr, 
and the magnitude V(£) — Vr is at most AV-. We will explicitly see that the two terms 
in the second part, which correspond to 5r and <ft contribution are of the same order and 
therefore the computation of 6r(£) in Eq. ( 4.12| ) is important. The last part represents all 
the high order terms of AV-/Vr- Notice that by the perturbation 5r, ^ m ax is no longer 
tt/H. However, the action change induced by (£ ma z — jj) is of the second order of the 
AV-/Vt and hence is included in ( 4.15| ). 

In conclusion, if we can get the analytic solution of (4.9) then the leading order of r 
is immediately given by the integral in ( 4.1 2| ) . It then follows that ( |4.13j ) and (4.14) give 
the leading order Euclidean action. Subtracting Se{4>+) given by([4.8|) from it, we get the 
factor B. 



5. The Four Scenarios of Tunneling in de Sitter Space 

In general, Eq. is still too difficult to solve in a way that the physics is transparent. 

However, for the triangle potential, we can solve ([Tl]) analytically. For the different choices 
of the five parameters Vt, Acp + , A</>_, X±, we find here are four different kinds of Euclidean 
solutions for <fi (shown in Figure |I]) : 

• Case (I) </>(£) reaches both <p + and </>_. This happens when a\/H 2 A(j) + > 1 and 
(\/H 2 )l(H 2 A(j) + /\^ > A(f>- where a ~ 0.4 and I(x) is a monotonic function shown 
in Figure ^. 

• Case (II) (/>(£) reaches cj) + but not That is, 0(^) reaches 0+ and (ftt-, where 4>t < 
(f> t - < 0_. This happens when both a\/H 2 A(t) + > 1 and (\/H 2 )l(H 2 A(f) + /X) < 
A0_ are satisfied. 

• Case (III). </>(£) reaches (ft- but not (ft + . That is, i^(^) reaches </>/+, where cftr > (ftf+ > 
(ft + . This case can happen only when the gradient towards the true vacuum is steeper 
than that towards the false vacuum. 
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• Case (IV). reaches neither <p + nor <^>_. That is, </>(£) reaches only 0j + and <§>%—. 
This case happens when a\/H 2 A<p + < 1 is satisfied. 

In the absence of gravity, the triangle potential just has two kinds of bounce solutions 
0, the cases (I) and (II). Here, including gravity, we have a new bounce solution, either 
case (III) or case (IV) , because here there is a cut-off £, m ax of the radius coordination £ due 
to the Hubble radius of de Sitter space. So in the finite range < £ < ir/H, the bounce 
solution may reach neither <j) + nor </>_. We will see that case (IV) is like HM tunneling, 
the case (I) is like thin-wall CDL tunneling, and case (II) is an intermediate case between 
the two limits. By adjusting the five parameters of the potential, we get a transition from 
HM tunneling to thin-wall CDL tunneling. 

For the triangular potential V((f>) shown in Figure ||, we have, for large Vt, 

4>" + 3H cot(#£W = ^r = ±A± (5.1) 

d(p 

whose general solution is given by 

HO = ai ± ^ + hf 2 (H0] (5.2) 

where f\(HO is the special solution and f% is the homogeneous solution. 

3 11 

fi(x) = - cot 2 x - esc 2 x lnfsinx), 

JiK ' 8 24 3 v h 

Hx) = -icsc 2 (|) + Iln(tan(|)) + I sec 2 (|). (5.3) 

Note that f^(HO = H/ sin 3 (ff£). The constants ai and b\ are to be determined by the 
boundary and matching conditions. We begin our discussion with case (IV) since it is new 
and its limiting case corresponds to HM tunneling. 

5.1 Case (IV) 

For case (IV), the bounce solution of ( |5.l|) contains two pieces, which correspond to the 
(/>_ side and the <j) + side, 

where (j>- > (ftt- > <Pr(0 > 4>T and < <pf + < 4>l{0 < ^T, whose general form can be 
obtained analytically. The boundary and matching conditions are 

^h(O) = o 

0r(£t) = 4>lHt) = 4>T (c . ,v 

tffettr) = </>' L (fr) 1 ' J 

^(TT/ff) = 

The first and last conditions require the specific combination fx + 2/2/3 in the solutions 
Q5.2| ). Imposing the boundary conditions at we have, 

MO =<t>T- ^{f(HO - f(H£ T )) 

MO = + |f (/(tt - #0 - /(vr - fffr)), (5.6) 
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<t>- 




Figure 4: The Euclidean bounce solution of </>(£) as a function of £ = [0, tt/H] in the four cases: 
(I) 4> reaches the true vacuum at 0_ inside the nucleation bubble and the false vacuum at cj>+ 
away from the outside of the bubble, (II) reaches </>+ but not (/>_ , (III) <fi reaches but not cf> + 
and (IV) <p reaches neither 0_ nor (f> + ; here the bounce starts at (f)(0) = cf>t- < (f>- and ends at 
4>(£.max = Tt/H) = <fif+ > cf>+. Case (I) and (II) can happen in the absence of gravity. Case (III) and 
(IV) can happen due to the presence of the de Sitter horizon. The thin-wall approximation is valid 
when the transition from 0_ to </>+ in case (I) is rapid. The HM formula is a good approximation 
when the variation of deviates little from <f>T in case (IV). 



where 



ax) = h{ X ) + 2 -h( X ) =yS asgc (i y - 8in ( 2c ° s (i ) 2 ) - 9 1 • (5 - ?) 



where we note that /(tt — x) = f\(x) — 2/2(x)/3. Matching the derivative at £r determines 

A_ f'(ir-H£ T ) 

c= a7 = Jvm (5 - 8) 

so all parameters are now fixed. 

To simplify the discussion, let us first consider the symmetric case A = A + = A_, so 
the analytical solution has a symmetry about £ = 7r/(2if), and the solution is particularly 
simple: £t = n/(2H). It is easy to check that </>(£) is a monotonic decreasing function. For 
this case (IV), we require that does not reach 4> + or in terms of the solution (|5.6|), 

/+ = <p L (ir/H) > </>+, (5.9) 

Introducing a useful dimensionless parameter 

^nkf + < 510 > 

the above condition means 

aA . 

" 7 = ww + < 1 (5 - n) 
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where 

a = /(tt/2) - /(O) = -1 + (A + !^ „ 0.398 (5.12) 

is a numerical constant. The similar condition 4>r(0) < <ft- is satisfied automatically since 
A</>_ > A0 + . In conclusion, if the parameters of the potential satisfy (5.11), then the 
bounce solution is case (IV). 

Before computing 5r and 5^(0), we estimate the magnitude of (/)(£,)■ By (| 



When A is small and H is large, 0(£) is confined in a small region which centers at 4>t- 
Because (ft = 4>t is the Hawking-Moss bounce solution, case (IV) tunneling is like a fluctu- 
ation around HM tunneling and will approach the Hawking-Moss solution when X/H 2 is 
small. 

Further, for the (shifted) Euclidean energy L, 

2 (5-14) 



V H 2 ' 1 J| H 2 

Hence L/F T ~ X 2 H- 2 /V T . By the condition Q and ( CT) , the ratio X 2 /(H 2 V T ) < 1 
and hence L/Vr <C 1 and the expansion ( pl|) and ( ETUI) works for case (IV). The only 
new issue is that for case (IV), the small expansion factor is X 2 /H 2 Vr not AV-/Vr, simply 
because the bounce solution ( |5.6| ) does not "feel" (ft-. The integral (4.12) is carried out 
analytically, and by ( pl| ), pl^ ) and ( |4~T5| ), we find that the two terms in the leading 
order correction have the same form, 

24vr 2 M 4 X 2 M® /A 4 M D 8 \ 
SeW = " C—^ + 0Ui . (5.15) 



where C = C\ + Ci ~ 65.49 is a positive constant. C\ ~ —196.47 corresponds to the 5r(£) 
contribution and C2 ~ 261.96 is from (ft((,)- So we see explicitly that the two terms in Eq. 
( 4.14 ) are comparable. That is, the back-reaction is important. The above result has a 
simple interpretation that we shall describe in Sec. 6. 

Since the Hawking-Moss scenario is a special limit in this case, let us treat this decay 
as a HM transition, with a correction to the HM formula that can be explicitly evaluated. 
So B is given by 

„ A ( 1 1\ X 2 Ml fX 4 M%\ 
B = S E ((ft) - S(cft + ) = 24vr 2 M p 4 ^— - —J - C-^ + o(-^J. (5.16) 

The first term is just the Hawking-Moss bounce Bhm ( |1.4| ). while the second term, correc- 
tion from this bounce solution, lowers B and therefore the actual tunneling rate is larger 
than that given by the HM formula. 

For Vr >V+, the first term in the bounce ( |5,15D is largely cancelled by S((ft+), so the 
corrections can be important. Rewriting the bounce formula as 

c (v+\ \ „ A JV+ 

l\ - 



B = B HM ( 1 - ^ ) + ■■■) = B HM ( 1 - 0.2317(a 7 ) I ^ ) + .... ) (5.17) 
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Figure 5: The function g(c). 



Since both 0:7 < 1 ( 5.11| ) and V+/Vr < 1, we see that the correction to the HM formula is 
at most 23%. 1 



The picture of case (IV) is consistent with the interpretation proposed in [12], [13] 



1 14]. HM tunneling should be interpreted as a quantum fluctuation up the potential barrier. 
Here we see that the field will thermally fluctuate from V+ at </> + up to V{4>l{ , k / H)) < Vr 
at (f>L(ir/H) > 0+ j tunnel to <Pr(0) < 0_ and then classically roll down to </>_. In the 
limit where 7AV4 /Vr — » 0, Vt is simply reached via the thermal fluctuations alone. Here, 
we see that the correction is typically not very big and the HM formula is quite good in 
general. On the other hand, as 7 — > 1/a = 2.51, 4>l{^/H) — ► + ; the outside of the bubble 
can reach the false vacuum within the horizon. When this happens, the HM formula is no 
longer accurate. 

Now let us consider the more general case where c = A_/A+ 7^ 1. The solution of £r 
should be determined by 

c( - Mr) - = /i(£t) - \m T ) (5.18) 

which can be solved numerically. For large c, the solution of £t is determined by series 
expansion, £t ~ ^/16/ (3c). The condition of the case (IV) is now 

a(c) 7 < 1, (5.19) 

where a(c) = (/i(Ct) - f/2(6r)) - (/i(tt) - |/2(tt)), where a(l) = a ~ 0.398. For large c, 
by the asymptotic form of £r, a(c) ~ y / C(c)/12. We can integrate the solution to get Sr 



1 Notice that the perturbation series (4.14), (4.14) and (4.15) start from the de-Sitter space Euclidean 



action, not the HM tunneling exponential factor Bum- So the correction, though large compared to Bhm, 
is still much smaller than the zero order, de-Sitter space Euclidean action. Hence the leading correction 



yields a very good approximation here. That is, the last term in Eq.(5.15) or (p.lq) is negligible 
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and B, 

E5 M ^l- 9 (c)( a (c) 7 )|j. (5.21) 

As before, both a(c)j < 1 and V + /Vt < 1. The function g(c) is plotted in Figure ^, When 
C is large, C(c) grows as and numerically the proportional coefficient is ~ 131. Hence, 
by the asymptotic form of a(c), 

g(c) -> 0.64, when c oo. (5.22) 

which means that the leading order correction is at most about 64%. We see that the HM 
formula is an over estimate of the value of B, by as much as a factor of three. 

5.2 Case (III) 

Here we shall simply show that this case, where </>(£) reaches (/>_ but not <j) + in the Euclidean 
solution, exists for some choice of the potential. Following the above solution ( |5.6| ), we have 

A<t> R = </> R (0) ~<h = ^2 (/(#£r) - /(0)) < A0_ 
A0 L = 0r - 0l(tt/F) = ^ (/(vr - Ffr) - /(ir)) < A0+ (5.23) 

Let AV R = V T - V((j) R (0)) < AV_ and AV L = V T - V(</> L (ir/H)) < AV+. Case (IV) 
corresponds to 

AV R < AF_, AV L < AV+ (5.24) 
while case (III) would be reached if 

AVr = AVL, AV L < AV + . (5.25) 

Using Eq.( |5.7| ), Eq.( |5.8| ) can be rewritten as 

c = A_ = 2cos 2 (^T/2)+cot 2 (gg T /2) 

A+ 2sin 2 (tf£ r /2) +tan 2 (H£ T /2) 1 ' ' 

Moving slightly away from the symmetric case, we have, from the above equation ( 5,26j ), 
for small deviations, 

HZ T ~n/2-±(c-l). (5.27) 

so we have 

AV R cA(()r ,9a-2(c-l) , 4 W 

Hence for c > 1, we have 

AVr > AV L (5.29) 
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For AV- > AVr and AV+ > AVl, we have still have the case (IV). Now we can consider 
potentials where 

AV R > AV- > AV+ > AV L (5.30) 

(decreasing A0_ while keeping other parameters fixed) which belongs to case (III). This 
demonstrates that case (III) exists for some types of potential. Here is an explicit example: 
when Acf) + = 0.005, A<£_ = 0.004, A_ = 1 x 10~ 7 , A+ = 0.5 x 10" 7 , V T = 0.77 x 10~ 6 , the 
bounce solution, which is plotted in Figure.^ does reach <p_ but not </>+. 

On the other hand, for c < 1 (that is, when the gradient towards the false vacuum 
is steeper than that towards the true vacuum), the Euclidean solution has AVr < AVl- 
This means the tunneling is from Vl = Vr — AVl going up to Vr = Vt — AVr. As noted 
earlier, this tunneling up phenomenon does not happen in the absence of gravity. So this 



possibility is a gravitational effect. This phenomenon has been studied in 15]. 

Note that this case which reaches 4>- but not </>+ does not exist for the symmetric 
(A = A+ = A_) case. In this case, c = 1 so that AVr = AVl and the above condition 



Q5.30 ) cannot be satisfied unless AVr = AV+ = AVl, which means that <fi does reach 



the false vacuum value. This corresponds to the cases we shall now turn to in the next 
subsection. 

5.3 Case (II) 

To see the existence of this case, we may restrict ourselves to the symmetric (A + = A_) 
case. Since A^>_ > Ac/> + we can have the situation where 

A^>^> A<P + (5.31) 



That is, the condition (|5.1l[) does not hold. Instead, we now have 



a 7 = a^A_ > !. (5.32) 

In this case, </>(£) reaches the false vacuum V+ at <f>+ outside the bubble within the 
horizon but the bubble inside never reaches the true vacuum V- at (/>_ at the moment of 
creation. In this case, the solution of Q5.1|) contains three pieces. 



4i(0*€E[(Ut] 

M£)ze (5.33) 



where <fo.(£) > 4>T and 4>2(0 < <t>T- The general solution for <f>\{£) and 4>2{£,), which satisfies 
0i(O) =0, is 

MO = ^(C2 + h(H0 + Af 2 (H0 ) . I .-,.:{ I i 
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where f\ and f 2 are given by Eq,( |5.3[) . The remaining boundary and matching conditions 
are, 

h{ir) = 4*2 (£r) = 4>t 

^i(fr) = <#>(&) (5-35) 
^ b 2 (Z + ) = 4 + , = 0. 

We have five parameters to determine, C\, C 2 , A, £t, and also five equation from ( 5.35j ). 
It is helpful to rewrite ( |5.35| ) as three equations only in A, £t, £,+ 



/i(tffr) - h(H£ + ) + A(f 2 (H£ T ) - f 2 m+)) = H 2 A<p + /\ = l/ 7 
fi(HZ T ) + (A/2 + l/3)&(HZ T ) = 
{f[(HS + )+Aft(HZ + )=0 



(5.36) 



Because the functions fi and f 2 contain only pure number coefficients, the dimensionless 
parameters A, H£t, H£+ are determined just by the combination 7 = A/ ' H 2 : A0+. So we 
may rewrite them as functions of I/7 which can be obtained numerically. Then we can 
insert A, H£t, back into ( 5.35| ), we will get all the parameters. 
For case (II), </>(£) does not reach <f>_, or 



or 



l(O)-0 t < A0_ 



hm T ) + \f2(Hi T ) - /i(0) - ~/ 2 (0) = J(l/ 7 ) < 



(5.37) 



(5.38) 



Notice that the left hand side is just a function of I/7 = H 2 A(ft+ / \, which we call I(l/A). 
So we need, 

'H 2 Afa 



_\_ 
H 2 



A 



< A0_ 



(5.39) 



Equation ( |5.39 ), together with ( f>. 32 ), is the condition for case (II). The function I(x) is a 
monotonic function which is plotted in Figure.0. Notice that I(x) > x, therefore, 



(f>(0) -<f>T > A0., 



(5.40) 



which means the quantum tunneling is from the false vacuum V+ to someplace lower than 
V+ but not exactly at V- . After the realization of quantum tunneling, V will continue to 
drop classically until reaching V-, the true vacuum. After getting the solution (ft, one uses 
the formulae (|4.12Q , ( 4.13| ) and ( 4.14 ) to get the factor B for the case (II). Notice here 
(ft does not reach </>_, so actually, the small expansion factor is not AV-/Vr but an even 
smaller factor AV + /Vt- 



5.4 Case (I) 

In the this case, (/>(£) reaches both q 

c*7 
A 



and <p_. Therefore the condition is 
A 



a 



H 2 



H 2 A(ft 
H 2 A(ft+ 
A 



> 1 



> A(ft_. 



(5.41) 
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Kl A) 




Figure 6: The function 1(1/ j). Note that 1/7 < a = 0.398. 



0(£) contains four pieces, 



HO 



4_ a €[<),£_] 

*> 2 (0 x e 

where 0i(£) > 0T and 2 (£) < 0T- The general solution for 0i(£) and 2 (£) is, 

A 



(5.42) 



MO = jp\ c i- hint) + A 1 f 2 (H0 



fc(0 = Jp[C2 + h(H0 + A 2 f 2 (H0 ) , 



(5.43) 



where /1 and f 2 are defined in (|5.3| ). The boundary conditions are, 



^i(e-) 

01 (£r) 

0i(fr) 
I 2 (£+) 



02 (£r) = 

2 (fr) 
0+) 02 (£h 







(5.44) 



So we have seven parameters to determine, C±, C 2 , A\, A 2 , £+, and also seven equa- 

tions from ( 5.44| ). Like case (II), we can reduce the number of both the parameters and 
the equations and then solve it. The solution looks like the thin- wall tunneling solution 
in Q. The reason is, by the condition (5.41), A<fi+ cannot be very small comparing with 
A0_, so the energy difference of the two vacua e = V+ — V_ cannot be large and the 
thin- wall approximation may apply. Again, after getting the solution of 0(£), we can get 
the Euclidean action and B by (|4.12| ), Q4.13D and ( jig) . 

It is clear that the thin- wall approximation belongs to case (I), since the inside of the 
bubble reaches 0_ while the outside reaches 0+ (see Figure ||) . The thin- wall approximation 
requires a rapid transition of from 0_ to 0+. This case has been analyzed in Ref. ||, 
1 16], 1 17]. For the sake of completeness, let us review the basic result following Q. In 



- 21 - 



the thin- wall approximation (e — > 0), we may divide the integration for the bounce B into 
three parts. Outside the bubble, <\> = 4> + and thus 



B„,,t = 0. 



In the wall, we have 



B wa u = 2ttVt 



(5.45) 



(5.46) 



where r is the bubble size and r is the tension of the wall which is decided by the barrier 
between the false and true vacua, 



Inside the bubble, 



Hence 



is a constant and Eq.(f2~i|) becomes 
d£ = dr{l-Kr 2 V/3)- 1/2 . 



S 



12vH 



E,in{ 



K JO 

Summing the three parts of B, we obtain 



rdf{\ - nV((p)r 2 /3) 1/2 . 



B = 2ir 2 r 3 T + 



i-ftl-^./S)* 



f ((1 - «rV + /3) 



3/2 



The coefficient B is stationary at r = R which satisfies 



R 2 ~ 9r 2 



k(V + + V-) k 2 t 2 



6 



+ 



16 ' 



where e = V+ — VL . So 



B tw = 2n 2 R 3 r + 



At: 2 



R 2 



'1-H 2 R 2 ) 3/2 -1 



Rlin-HiR 2 ) 3 ' 2 -! 



(5.47) 

(5.48) 
(5.49) 

(5.50) 

(5.51) 
(5.52) 



where H^. 1 = R± = (kV±/3) -1 / 2 and the subscript "tw" means the thin- wall approxi- 
mation. Note that the last term is proportional to k 2 and so is very small most of the 
time. According to Eq. (|5.51| ), we can easily check that the bubble radius at the moment 
of materialization is not larger than the event horizon R + of the de Sitter space in false 
vacuum. This is reasonable; otherwise the bubble cannot be generated causally. 

For the special case with R <C R+ , the bubble size is much smaller than the curvature 
radius of the background and gravity does not play a big role. In this limit B becomes 



B = 2ir 2 r 3 T - — r 4 e. 

2 



Here, B is stationary at 



r = Rn 



3r 



(5.53) 
(5.54) 
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0.05 0.10 0.15 0.20 

e 



Figure 7: The contour plot of (Bt w — B)/B as the function of e and Vt- The horizontal axis is e 
and the vertical axis is Vt- The contour line labels the value (B tw — B)/B and the darker regions 
have larger differences. 



SO 

27tt 2 r 4 , . 

Btw ~ -= r ■ 5.55 

2 

In the other limit of R ~ i? + which corresponds to V+ ~ VI = V 3> Vs = 2e 2 /3KT 2 + 
3kt 2 /8. This happens at the high energy scale in the landscape, and the bubble radius is 
given by R ~ w3/ 'kV . Now B is dominated by the first term in Eq.( ^.5Q ), namely 

B tw ~ 6V3vr 2 r(KV)- 3/2 = (5.56) 

In the units where M p = 1, or equivalently k = 1, the tension of the bubble satisfies r<l. 
In Planck region (V ~ 1), B <C 1 and T ~ 1. At low energy scale, the background curvature 
radius is quite large and the bubble size is relatively small, and then the tunneling rate is 
insensitive to the vacuum energy of the false vacuum. 

The familiar CDL thin- wall approximation occurs in case (I), however, a solution 
of the case (I) is not necessarily the "thin- wall" solution, i.e., the tunneling factor Bt w 
determined by the thin-wall approximation may not be very accurate. To illustrate this 
point we compare the values of Bf W with our computation for case (I). 

We fix G = 1, A = 1 x 1CT 7 , 4>t = 0, (f>+ = —0.005 and vary the variables <p- and Vt- 
We define 

« v.-v- 0-1^1 (! . 57) 



AV+ V T - V 
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which is always positive and we may expect that when e is close to 0, the thin-wall ap- 
proximation is accurate. We reserve B for our result for the tunneling exponential factor 
from Eq.( 4.12| ), Eq.( 4.13| ) and Eq.(4.14) while Bt w denotes the counterpart of the thin-wall 



approximation formula ( |5.52 ) given in Ref.Q. We plot the relative deviation of the two 
methods, (Bf W — B)/B, as a function of r and Vr in the contours of Figure [/]. 

In the contour figure we draw the range 0.01 < e < 0.2 and 5 x 10~ 8 < Vr < 7.5 x 10~ 7 . 
The reason for the Vr range choice is that if Vr is too small then our expansion is not good 
and if Vr is too large, case (I) will turn into case (II). From the figure we can see 

• The contour lines are roughly vertical and the relative difference (Bf W — B)/B in- 
creases from the left to right. It means that when e is small the thin- wall approxima- 
tion coincides with our computation and so verifies that the thin-wall appoximation 
is good when e is small. 

• However, the contour lines are not completely vertical but tilt from the top left 
to the bottom right. It means when e is fixed and Vr is increasing, the difference 
is getting larger and larger. An immediate explanation is: in this situation the 
solution is moving toward case (II) and finally when tunneling becomes case (II) 
the difference is very large, since the solution does not reach 0_ in contrast to the 
thin-wall approximation. Or by the Euclidean equation of motion, if Vr is large, H, 
which serves as a damping term, is also large and slows the evolution of the Euclidean 
solution. So the wall is not "thin" as before. 

Although unlike case (IV), it is hard to get the analytic form of B for case (I) because of 
the transcendental Eq. fl5.44 ). However, it is helpful to look at the numerical computation. 



For example, for a potential with small e = V+ — V_, like A = 1 x 10 7 , Vr = 0.2 x 10 6 , 
A0+ = 0.005 and A0_ = 0.0055, the zero-order of the S E , by Eq.(PH), is S E ,o = 



-1.87500000 x 10 b . Eq. (|4.14|) gives the leading order correction S E ,i = S E ,i,6r + S E ,i,, 
where Se^St = 4787.05 is the first term in ( 4.14 ) from 5r contribution and S E ,i,<j> 



—7691.97 is the second term induced by (p. Therefore, we explicitly see that the leading- 
order corrections S E ,i,Sr and Se,i,4> are much smaller than the zero order S E ,o and the 
perturbation series for S E is valid. To get the factor B, we needs to compute S(<fi+) which 
is -1.8796992 x 10 6 , so 

B = S E - S((p + ) ~ Se,o + S E ,i,s r + S E,i,<t> ~ S((j>+) 
= 1794.32 

which is quite close to the thin- wall approximation result Bt w = 1848.76 by Eq.( |5.50| ). 

We see that it is crucial to include the back-reaction effect due to 5r here. In fact, if we 
drop the contribution from 5r in Eq.( 4.14| ), then f?((5r-excluded) ~ S Et o + Sp;^^ — S(4> + ) ~ 



—2992, which is negative and has no physical meaning. The importance of the Sr term in 
Eq.( [4,14 ) cannot be over-emphasized. 
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6. Discussion 



Let us comment on the physical meaning of the cases other than case (I). For case (II), 
after the realization of quantum tunneling, the field will evolve by the classical equation. 
Coleman [|| points out that, instead of going back to solve the Minkowski (Lorentz) equa- 
tion, we can get the classical solution directly from the Euclidean solution by analytical 
continuation. For ||, we can verify that, after the quantum tunneling, the field will con- 
tinue to roll down until reaching the true vacuum. This also happens in de Sitter space. 
After the creation of the nucleation bubble with radius R at time t = 0, the evolution of the 
surface of the bubble is described by the Lorentzian action, where £ 2 = R 2 — > r 2 — t 2 . The 
center of the bubble starts at £ = and as t increases, £ 2 = — t 2 becomes negative at the 
center. By analytic continuation, starting with (ft{r = 0, t = 0) < (ft—, <p(0,t) falls towards 
(ft-. It may overshoot and oscillates about the true vacuum at (ft-. With Hubble (or any 
additional) damping, it will eventually settle at the true vacuum as the nucleation bubble 
continues to grow. In case (III) and case (IV), (ft at the bubble outside is also expected 
to follow classical motion and roll down towards the false vacuum at (ft + (while in a GH 
thermal bath). However, since the bubble is growing rapidly at the same time, the true 
vacuum may be reached before (ft has time to roll to the false vacuum. 

Let us show the different features of the tunneling by an explicit example. The type 
of tunneling is determined by the parameters of the potential, Vr, A+, A_, 6<ft+ and 8(ft-. 
First, we want to see the dependence of the exponent B on H, for large H. We set G = 1 
and fix <ft T = 0, A(ft- = 0.007, A(ft+ = 0.005, A± = 1 x 10 -7 so the shape of the potential 
is fixed. We can move the potential up and down by varying Vr (hence H) to see the 
dependence of B on H. 

• The condition for case (IV), 07 < 1 reads Vr > 9.49476 x 10~ 7 . For example, we can 
plot the solution of cft{r) for Vr = 2 x 10~ 6 in the last picture of Figure ^| and it is 
clear that neither <ft + nor (ft- is reached. The exponential factor B can be obtained 
by Eq,( |5.16 ), B ~ 41.7301. The HM tunneling exponential factor for this potential 



is Bhm ~ 46.8867 hence the case (V) tunneling rate is close to the HM tunneling 
but faster than it. The relative difference of B is about 11%. 

• When aj > 1, we need to consider the condition 

-^1 — ^± <A0- (6.1) 



H 2 V A 

which reads 9.49476 x 10~ 7 >V T > 6.46586 x 10~ 7 . If it is satisfied, the tunneling 
is of case (II). For example, we plot the solution for Vr = 9.0 x 10~ 7 in Figure |I| and 
find that (ft reaches but not (ft-. 

If this condition is also violated, i.e., Vr < 6.46586 x 10~ 7 , the tunneling is case (I). 
2 For example, the solution of Vr = 3.0 x 10~ 7 is plotted in Figure [| and (ft reaches 
both (ft + and (ft-. Eq.( 4.14| ) gives B ~ 928.714 while the thin-wall approximation 



^Notice that our approximation works for Vr S> AV_, which reads Vr ~S> 7 X 10 
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Eq. (|5,52| ) gives B tw ~ 1019.82. The relative difference is about 9.8% due to the finite 
difference of V+ and V-. (So the bubble wall is not "thin" enough for this potential.) 

As H is decreased from a large value, the tunneling type goes from case (IV) to case 
(II) and finally to case (I). The corresponding exponent B for different cases are plotted 
in Figure |as a function of H. It is straightforward to see the dependence on the shape of 
potential, say, by varying A = \± and fixing A</>_ = 0.007, = 0.005, Vr = 8 x 10~ 7 , 
we can also see the transition among the three cases. Notice the combination xl(l/x) is 
monotonically increasing function of x, so the conditions can be solved easily in A, 



• For A < 8.4257 x 10~ 8 we have case (IV). 

• For 8.4257 x 10~ 8 < A < 1.23727 x 10" 7 we have case (II). 

• for A > 1.23727 x 10~ 7 we have case (I). 



Therefore case (IV), which is close to HM tunneling, happens for a "flatter" barrier (small 
A) while case (I), of which thin- wall tunneling is a special case, happens for a "sharper" 
barrier (large A). 

We may write the factor B ( |1.2[ ) as composed of two terms, 

B = [S E (<f>f+) ~ S E (ct>+)] + [SE^bounce) ~ S E ^f+)\ (6.2) 



where (j)T > 0/+ > (f>+- For cases (I) and (II), we have 0/+ = </>+, so the first term vanishes. 
For cases (III) and (IV), (pf + > + , so both terms contribute. In the limit cf>f + — > <j>T, 
</>(£) — > <f>T, so the second term vanishes and the resulting formula reduces to the HM 



limit. The picture is consistent with the interpretation in [12, 13, 14 1. However, as we 
have pointed out, corrections and back-reaction can introduce large corrections to the HM 
formula. 

Brown and Weinberg || showed that the CDL tunneling rate can be derived by treating 
the field theory on a static patch of de Sitter space as a thermal system. In this thermal 
system, tunneling does not need to occur from the bottom of the false potential well. 
Instead, tunneling proceeds by a combination of thermal excitation part way up the barrier 
followed by quantum tunneling through the barrier. The tunneling rate is a thermal average 
of the energy-dependent quantum tunneling rates 1£], 

' dEe -{E-E{4> + ))/T e -J{E) (6 3) 

where E is the energy of the field configuration. Here the quantum tunneling rate is given 
by WKB approximation J{E) = 2 f£*~ d(j)y/2(y{<p) — E) where the integral is along a 
Euclidean path from one classical turning point <pf + to the other classical turning point 
<frt-- This integral is dominated by the energy E* that maximizes the integrand. Using the 
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saddle-point approximation it follows that 



i _ _ r l - , , ! 

T ~ 



h+ y/2(V(<t>)-E) 
r<h- l 

= 2 / d^-j— (6.4) 

= 2Ar (6.5) 

since the integral is along a solution to the Euclidean equations of motion. Fixing T = Tjj 
determines E* (and so 0/+ and <f>t-)- This yields J{E 1f ) = SE{4>bounce) — ££(</>/+). This 
result also implies that the tunneling takes a Euclidean time At = l/2Tfj and 

i 

* E(<P)dT ~ E/T H (6.6) 

That is, E^/Th = Se{4>}+) an d E(4> + )/Th = Se(4>+)- Using this result in the saddle-point 
approximation to ( |6.3[ ), one finds 

r ~ e -(s E {<t>f + )-s E (4>+)) e -{s E {4> bounce )-s E {4>f + )) (6.7) 



which reproduces the standard CDL tunneling rate. So our result ( p.2|) agrees with this 
result of Ref.[||: the first term in ( |6.2j ) corresponds to a thermal (i.e., Gibbons-Hawking 
temperature) fluctuation from the false vacuum (at </>_)_) part way up the barrier to 4>f + and 
the second term in ( |6.2| ) corresponds to a quantum tunneling. Note that the derivation of 
Ref.|| assumes a fixed de Sitter background., while the back-reaction (which can be large) 
is included in our derivation. In this sense, our result is a big improvement. In Sec. 5.1, 
we calculated the S ' Ei^bounce) and determined for a given triangular potential. This 
allows us to compare quantitatively the contribution from the GH thermal effect to the 
contribution from the quantum effect in such a tunneling. 



7. Thermal Tunneling 

At finite temperature, the transition from V + to V_ can also occur via thermal fluctuations 



1 2 if] . So there is another contribution due entirely to the GH temperature fluctuation. The 
dominant thermal tunneling process is the formation of 0(3) symmetric bubbles of V- that 
minimize the change in entropy AS. Let us evaluate the rate of this process and compare 
it to the above tunneling process. The thermal tunneling rate is given by 

e AS = e AF/T {71) 

where AF is the change in the free energy. We can think of the free energy as a three- 
dimensional action 

AF = S 3 = [ d 3 x(l{d^) 2 + V(cj>)) (7.2) 
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for a scalar field with a standard kinetic term. In this case the equations of motion are 

2r' dV 

<p + — = — 7.3 

r acp 

We restrict our attention to the symmetric triangle potential as before. In the large 
Vt limit considered above, the general solution to (|7.4p is 



4> = ^{A+^-cot{Hi) + Bcot{Hi)) (7.5) 

Like the Euclidean case discussed above, there are also four types of thermal tunneling 
solutions: 

• Case I. reaches both <p + and 0_. 

• Case II. </>(£) reaches <p + but not 

• Case III. </>(£) reaches 4>- but not <p + . 

• Case IV. reaches neither + nor 0_. 

We will begin our discussion with case IV as it is the simplest. 
7.1 Case IV 

In this case, we need to solve ( |7.3| ) with the boundary conditions ( |5.5| ). The solution to 
these equations is 



i(0 


= 4>T + 


2(0 


= 4>t- 




7T 


6r 




~ 2H 



(7.6) 



This solution is valid if 



A^cot^) 

(7.8) 

< 1. (7.9) 



2A 

Comparing ( |7.9| ) to ( 5.1 1| ) we see that the existence of a case (IV) Euclidean tunneling 
solution implies the existence of a case IV thermal tunneling solution and that some poten- 
tials will admit a case IV thermal tunneling solution and a case (II) Euclidean tunneling 
solution. 

Evaluating ( |7.2| ) using Q7,8| ) we find that, for a specific background temperature, 

8irAV+ „/A 2 M f 6 

4 
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To compare to the HM bounce (|l.4| ) we find the ratio at leading order 

~ (7 - n) 

Based on the discussion in the last section, this result is expected and consistent with 
1 18]. One should compare this thermal tunneling to the stochastic tunneling process of 



7.2 Case II 

Case III tunneling only occurs when A+ 7^ A_ as in the Euclidean case. Since it does not 
show any special new feature, we shall skip this case. 



MO = 4l( A 2-^- cot(i?£) + B cot(fT£) ) (7.13) 



In case II, we need to solve ( [73| ) with the boundary conditions ( 5.35 ). Imposing the 
condition </4(0) = we find 

M0 = ^(A 1 + ^cot(H0) (7.12) 

_A_ 
IP ^ 

(7.14) 

As in the Euclidean case we can eliminate A\ and A2 by rewriting the boundary conditions 
( 5.35 ) as 

cot(tffr) + ^± cot 2 (#£+) + 5(cot(^ T ) - cot(^+)) = (7-15) 
cot(F£ T ) - H^ T csc 2 (i/£ T ) - B csc 2 {H^ T ) = (7.16) 
\cot{Hi + ) - ^±csc 2 (^+) - Bc S c 2 (H£ + ) = (7.17) 

which we can solve numerically for B, and as functions of — ^ + . We can then use 
these solutions in the original equations ( |5,35 ) to determine A\ and A2. 

For our solution to be consistent, we need (f>i(0) — <pT < A<^_ or equivalently 



, H 2 A<f> + \ 1 Hfr ff 2 A0_ 
J s 3 ( ^ 1 = 2 2~ cot (^ T ) < a ' ' ' 



Numerically we find that Is 3 { H ^ + ) > ~ ^ + so case ^ thermal tunneling is always from 
to some Vp where V+ > Vf > V- as we saw above in case (II) Euclidean tunneling. 

7.3 Case I 

When ( [79| ) and ( [7.18 ) are both violated we have case I thermal tunneling. Here we need to 
solve ( |7.3| ) with the boundary conditions (5.44) numerically if the thin- wall approximation is 
not valid. In general the solutions have the same qualitative features as in case I Euclidean 
tunneling, but the Euclidean bounce B is smaller than S^/Th- 

Thin-wall CDL tunneling is a special case of I. The thin-wall approximation greatly 
simplifies calculations so for the remainder of this section we consider a general potential 
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that satisfies the thin-wall conditions. The difference between the CDL and thermal tun- 
neling rates has to do with R 4 space in CDL versus i? 3 x S 1 space in the thermal case, 
where S 1 corresponds to = 1/T = 2n/H. In the thin- wall approximation we can divide 
the integration in |7.2j into three parts. For a bubble of radius r we have S 3:OU t = and 
Sz,waii = 4:irr 2 T. Using Eq.( |5.48 ), the integral inside the bubble becomes 

SsmW = ~— f dH 1 ~ ^(</>)r 2 /3) 1/2 - (7.19) 



« Jo 

The free energy is given by 

S 3 = 4vrr 2 r+^ r (l - r 2 /R 2 + ) 1/2 + R + arcsin (j^j - r (l - r 2 /R 2 _ ) 1/2 - R arcsin 

(7.20) 

Maximizing S3, we find R given by 

1 1 (H 2 +Hl) K 2 T 2 

2 + ^ (7 - 21) 

(3) 

where R ' = 2r/e. In the R«i? + case; this yields 

5 3 = 47ri? 2 r - 47ri? 3 e/3 (7.22) 

(3) 

For critical size radius i? = R , we have 

Sg/T = 327T 2 r 3 /3e 2 F (7.23) 
To see which path dominates, we compare this to B = Se, 



B 27 R { n 4) 



(7.24) 



(syr) 64 r+ 

Since R <C R+, we see that the Euclidean S bounce B is smaller so it dominates. 
For large, R<R+, 

S 3 /T = 8tt 2 t/H 3 (7.25) 
so we see that B ( |5.56 ) is smaller than S 3 /T by a factor of four, 

B 1 



(S3/T) 4 



(7.26) 



It is easy to understand the origin of this factor of four. In CDL, we have 5* 3 with size 2-7T 2 
versus S 2 x S 1 with size 4-7T x 2ir = 8tt 2 in the thermal case. This just shows that the 0(4) 
symmetry lowers the Euclidean action B and so yields the correct answer. 
Note that this is true only for B>1. Recall that 

r = A e- B + A T e~ S3/T + ... ~ A x e~ B + A T e~ AB (7.27) 



- 30 - 



If this condition is not satisfied, say when B > 1, then this factor of four difference may 



be overcome by the difference in the prefactors Aq and At- Callan and Coleman [2C] show 
that, for bounce B, 

B\ 2 / det(-^ + V^ + ) \ V2 

2tt) V det'(-5 2 + V"{(j>) ) { ' 

where the prime on the determinant implies that the four zero modes are removed, yielding 
the first factor. This prefactor is difficult to evaluate in general, so here we shall use 
dimensional arguments to find an order of magnitude estimate. 

For B < 1, the above formula for CDL tunneling actually breaks down. Linde j2l| 
argued that 

3/2 



r ~ r4 G^) e ~ sz,T (7 - 29) 

The difference is that here, the temperature T is the Gibbons-Hawking temperature Th, 
which strictly speaking is a gravitational, not thermal, effect. 

This means that the potential V should now include the finite temperature effect: 
V((f)i) — ► V((pi,Tn)- The effective potential in de Sitter space has been calculated for some 
simple potentials in 0, M. As expected based on the GH temperature interpretation, the 
effective potential for scalar electrodynamics in de Sitter space calculated in |?J shows the 
same behavior as one varies the inverse radius of de Sitter space as the effective potential 
in Minkowski space as one varies the temperature. For a phenomenological potential (as in 
a slow-roll inflationary scenario), one may assume this finite temperature effect (and other 
quantum effects) is already built into the potential. However, in string theory applica- 
tions to cosmology, where the effective potential can be calculated given a specific model, 
such finite GH temperature effects should be included. Here let us consider some specific 
examples. 

8. Tunneling in the Cosmic Landscape 

We consider a type IIB compactification with the moduli stablized by a combination of 



fluxes and a nonperturbative superpotential as in [22]. These compactifications typically 
have a large number of axions (pi corresponding to integrals of the four-form potential 
over each of the independent four-cycles. Independent nonperturbative effects in the flux- 
induced superpotential give rise to a periodic potential 

V{cpi) = M 4 e~ s ^ ( 1 - cos ( ^ ) ) + V(p) (8.1) 



Ji 

where fi is the decay constant, M is a natural mass scale (say, the string scale), and 
the ith instanton has action S\ nst . V(p) contains the potential coming from the moduli 
and includes D-term contributions. The shift symmetries associated with these axions are 
broken by nonperturbative effects. Including the effect of a single Euclidean D3-brane for 
each independent four-cycle the axion potential becomes f23[ 

Vfa) = Vofa) + £ a t cos (|) + Y, A? ™ S (| - ^f) (8.2) 

i ij ^ 



- 31 - 



{Log 10 ^,Log 10 fi} 




Figure 8: Comparison of the Euclidean action and the free energy divided by the GH temperature. 
log 10 B (the black solid curve) and log 10 (5 , 3/Tf/) (the red dashed curve) plotted against log 10 V- 
with t = 10~ 10 and e = 1CP 22 (in reduced Planck units) constant. In the thin-wall approximation, 
CDL tunneling has an exponentially faster rate than thermal tunneling. We see that B varies by 
25 orders of magnitude. Tunneling is exponentially enhanced when the wavefunction is higher up 
(i.e., larger Vt) in the landscape. 



where Vo(fa) is a smooth function of fa due to D-terms. Here, typical a, are exponentially 
small compared to the string scale. It is easy to estimate t\ along the fa direction. Let 
«i = &i + ^j^Pij, so the height of the potential barrier in the fa direction is 2«j. Then, 
using Eq. fl5.47l ), we have 



Z<f>+ 




'2(2a i )cos(^) = - A /^/ i (8.3) 



and a crude estimate of the prefactor gives 

(I) 



r ^ ( \V'^t)\ 2 e- B (8.4) 



where 



and \V"(fa)\ = a/f 2 , yielding 



B = 2it 2 t/H 3 (8.5) 



vr 4 a 



i=i jt 

Tunneling can be fast if the axion decay constant / is sufficiently small. Cosmological 
upper bounds on the amount of axionic dark matter constrain / to be below 10 12 GeV, and 
astrophysical measurements of the cooling of red giants constrain / to be at least 10 9 GeV 



1 24]. For example if we take fa = 10 10 GeV within this observationally preferred region 
and choose a, = (2 * KT 7 M P ) 4 , then T ~ 1 for H = 1CT 7 M p and N ~ 10. For smaller / 
the tunneling rate T can be order one for even smaller H. For model-independent axions 
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in heterotic string theory / is generically between 1.1 * 10 16 GeV and the reduced Planck 
mass 



i \ I " I v\ (I'll .1 . .1. -i- 

\2§. If we choose ft > 1.1 * 10 16 GeV for N < 100 there is no choice of cti = a and 
H in the regime where effective field theory is valid (a <C Mp and H < M p ) so that T ~ 1. 
For relatively large H or T, we have, instead 



(^€\\V"(4>u)?e- S ^ T (8.7) 



where fat is the value of fa at the top of the barrier. 

So it is reasonable to expect that the wavefunction of the universe tends to spread 
along some of the axionic directions. This what we expect for the QCD vacuum; that is, 
the wavefunction is a Bloch wave with angle Oqcd- As we go up in the cosmic landscape 
(say, turning on D-terms), tunneling will be faster and so the wavefunction will be Bloch 
wave-like in more of the axionic directions. If we start with the wavefeunction localized as 
a classically stable vacuum site, it will take time for the wavefunction to spread. This time 
will be shorter when we are higher up in the landscape and when there are more axionic 
directions present. 

However, here V(fa) should be replaced by V(fa,T). Instead of calculating this, we can 
make another estimate. In quantum mechanics, tunneling through a barrier becomes less 
suppressed when the energy of the particle increases. This happens when the particle is in a 
thermal bath with rising temperature. In quantum field theory, thermal effects typically lift 
the potential in a way such that the tunneling becomes faster, until the barrier disappears 
(i.e., the tunneling probability approaches unity). 

For A0 4 theory, with 

V(faT = 0) = -^ + ^ 
we have, for high temperature T, 

V{faT H ) = \(^-rn*y + ^ + ... (8.8) 

so the effective mass term is no longer tachyonic for T > T c , where the critical temperature 
is given by T 2 ~ 24m 2 / A. 

For any direction in the moduli space, if 

24m 2 

T 2 H >— -S-, i = l,2,....,d (8.9) 

Expanding the above potential about a maximum (top of a barrier), we have 

a 2 a 
2Z 2 + 4!/4 



^ = -^ + 7774^ + 



we see that, if in any direction, 

24/ 2 < Tjj (8.10) 

then there is no barrier in that direction so the wavefunction is coherent along that di- 
rection. Even if this condition is not satisfied, we see that the finite Tjj effect will enable 
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to wavefunction to be coherent along some direction quickly, approaching a wavefunction 
similar to a Bloch wave. 

In fact, compactification and moduli stabilization usually introduces a term like H 2 (p 2 , 
which will tend to remove a lot of barriers. This term arises from Gauss' Law. Actually, 
this is the term that causes the rj problem for slow-roll inflation. 

In a single dimension at low scales these effects are probably negligible without a large 
degree of fine-tuning. For example consider a de Sitter vacuum in flux compactification 
in string theory with all moduli fixed. In the KKLT model |2^] , a highly warped type 
IIB compactification with nontrivial fluxes is stabilized using non-perturbative effects from 
Euclidean D-branes and gaugino condensation. The resulting anti-de Sitter minimum for 
the imaginary part of the volume modulus a is uplifted to a de Sitter minimum by adding 
a small number of D3 branes. The presence of the D3 branes induces a term D/a 3 where 
D depends on the warp factor and the number of D3 branes. The potential is: 

V= 2 a 2 ( s aaAe ~ aa + Wo + ) + ^ (8 ' n) 

In this model if one tunes the parameters to have the de Sitter minimum be the cosmological 
constant observed today the tunneling rate is exponentially long. 

The Gibbons-Hawking temperature effects are unimportant in the KKLT model be- 
cause in the regime where the calculation is under control (a 3> 1), the de Sitter vacuum 
is necessarily exponentially suppressed compared to the Planck scale. Increasing D both 
increases the false vacuum and decreases the barrier. The de Sitter vacuum disappears 
classically if the D-term uplifting is too large. Since the true vacuum is always Minkowski 
in this model HM tunneling dominates if the de Sitter vacuum is large. The HM bounce is 
B HM = 24tt 2 M*AV + /V£. Since M*/V T > 1 in this model, and AV + /V T <C 1 only when 
the barrier is about to disappear classically, the GH temperature corrections are small in 
this case. Numerical calculations approximating the maximum with a quartic potential 
and computing the finite temperature corrections using the GH temperature also show 
that the region of parameter space in which the de Sitter minimum is present classically 
but unstable quantum mechanically is extremely small. 

9. Summary and Remarks 

In this paper, we discuss the relation between CDL thin-wall tunneling and HM tunneling. 
The picture is in agreement with the qualitative understanding one has already, but the 
details allow us to obtain a quantitative understanding of the validity of each approximation 
and the sizes of the corrections. This result does allow us to say something about the stringy 
cosmic landscape. Here are a few comments that may be relevant, in no particular order: 
• When we are high up in the cosmic landscape, gravitational effects exponentially 
enhance the tunneling rate. Qualitatively, one may interpret this as a GH temperature 
effect, but the actually enhancement is exponentially larger than a naive estimate based on 
the usual temperature effect, due to the enhanced 0(4) symmetry of the Euclidean action 
(versus the SO(3) x U(l) symmetry in the finite temperature case). 
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• Above the GUT scale, quantum gravity effects become important and we have noth- 
ing to say in this case. Around or slightly below the GUT scale, the tunneling can easily 
be enhanced by so much that there is simply no exponential suppression at all and the 
semi-classical formula simply breaks down. This is in the single field case, with only one 
tunneling direction. In the multi-field (d > 1) case, as we expect to be the situation in the 
cosmic landscape, there are typically many tunneling directions for any false vacuum, so 
one expects fast tunneling to take place. 

• As is well known (see Appendix A in Ref. for a generic spherical potential in 
d > 3 spatial dimensions, the extra dimensions acts as an angular momentum like repulsive 
potential (with angular momentum I = {d — 3)/2). So it is harder to form a bound state 
in higher dimensions, as is the case in the cosmic landscape. Instead of having a bound 
wavefunction that has to tunnel out, we may have a resonance like situation where the 
decay time can be very short. Even if the wavefunction is trapped, it is more likely to be 
weakly trapped, so the wavefunction has a long tail outside the classically allowed region, 
rendering less suppressed tunneling. 

• It is again well known that the finite temperature effect on the effective potential 
tends to decrease the barrier, allowing faster tunneling. In fact, compactification and 
moduli stabilization usually introduces a term like H 2 (j) 2 , which will tend to remove a lot 
of barriers. This term, a finite volume effect, arises from the equivalent of the Gauss's 
Law. (Actually, this is the term that causes the rj problem for slow-roll inflation |2f|) Here, 
a similar term arises due to the finite temperature effect on the potential. One may re- 
interpret this effect as a finite volume effect, due to the presence of the Gibbons-Hawking 
horizon 1/H. 

• Even if the barriers are still present (that is, not completely lifted), we see that 
the finite Tjj effect will enable the wavefunction to be coherent along some direction more 
quickly. This is particularly likely along the axionic directions since the effective potential 
is periodic or close to being periodic and the heights of the barriers can be very low. In this 
situation, the wavefunction approaches a Bloch wave. Presumably, the QCD 6 vacuum we 
live in today is described by such a Bloch wave. 

• Since there are many vacua in the stringy landscape, resonance tunneling effects 
should play an important role in the landscape |2^] . This also enhances the tunneling 
among the vacua in the landscape. Such enhancement can be substantial. In particular, 
its role in the development of the Bloch wave is well known. 

• With all these observations, one suspects that the wavefunction of the universe in 
the cosmic landscape may be quite mobile; that is, it does not stay at any particular 
vacuum site long enough to allow eternal inflation. However, without a detailed knowledge 
of the structure of the landscape (vacuum sites and their nearby neighbors as well as the 
barriers between them), it is difficult to make definitive statements. As we just pointed out, 
the wavefunction most likely is spread out along the periodic directions (i.e., the axionic 
directions) as Bloch waves; so what happens to the aperiodic directions? If we treat such 
directions randomly, i.e., as a random potential, one may borrow the insight obtained in 
condensed matter physics to argue that the wavefunction should be fully mobile along these 
directions as well p7| . This implies that eternal inflation in the landscape is very unlikely. 
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This argument is based on the large dimensionality of the landscape and does not use 
the gravitational effects discussed in this paper; so it should apply all the way to very low 
vacuum energy density sites. 
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